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In this paper we utilize ^"-function regularization techniques in order to compute the Casimir force 
for massless scalar fields subject to Dirichlet and Neumann boundary conditions in the setting of 
the conical piston. The piston geometry is obtained by dividing the bounded generalized cone into 
two regions separated by its cross section positioned at a with a e (0, b) with b > 0. We obtain 

(N : 

expressions for the Casimir force that are valid in any dimension for both Dirichlet and Neumann 
boundary conditions in terms of the spectral ^"-function of the piston. As a particular case, we specify 
the piston to be a <f -dimensional sphere and present explicit results for d — 2, 3, 4, 5. 

I. INTRODUCTION 



The Casimir effect is one of the most important macroscopic manifestations of the zero point energy of 



X 



quantized fields under the influence of external conditions _l_l 



3811 or in spaces with non-trivial topology. In 



CM 

recent years, a vast amount of literature has been produced on the Casimir effect, which was first predicted 

aw y p jinn 

in the seminal paper [13], especially for its relevance in nanoscale physics H 10. 1 1 lL 13811 . Due to its nature, 



calculations of the vacuum energy lead to divergencies which need to be regularized and subsequently 
renormalized. Several regularization methods exist, amongst the most important ones are frequency cutoff, 



point splitting and zeta function regularization 



22, 



23 



38]. For many configurations, these 



techniques yield the same finite renormalized result, however the way divergencies are removed is different 
in each scheme. The non-uniqueness of the removal procedure raises the question, which of them is the 
physically best motivated one. Technical and interpretational problems of this nature can actually be avoided 
if one considers the Casimir effect between separate objects. In this case, the divergent part of the energy 
(for massless fields) depends on the heat kernel coefficient ao/2 related to the geometry of the objects. 
These coefficients, in turn, do not depend on the distance between the bodies and, hence, the Casimir force 
between them is free of divergencies Jioll . Belonging to the class of configurations for which the Casimir 
force has been unambiguously evaluated are pistons of certain types. 



These piston configurations, introduced in 11411 . have become increasingly important because of this 



fact. A large variety of piston configurations and boundary conditions have been studied throughout the lit- 
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erature, both at zero and at finite temperature. Rectan gula r Casimir pistons with different types of boundary 



conditions have been considered, for instance, in (21 
have been analyzed in [4J, 



32 



33 



37, 



simir pie 



Cylindrical Casimir pistons, instead, 



44I1 . Higher dimensional Casimir pistons with arbitrary cross sections have 



been studied in the setting of Kaluza- Klein models in [2' 



temperatures has been obtained, e.g. 



,inQ 



430. 



9. 



3511 . A generalization of these models for finite 



In the literature, Casimir piston configurations have the general geometric structure of a direct product 
I xN x M, where / is a closed interval of the real line, N is either a compact or non-compact manifold, and 
M is a compact manifold usually describing the extra-dimensions. Pistons of a different geometric nature, 
where the local line element contains a warp factor, have been recently considered in the framework of five 



26, 



4lh . Moreover, vacuum polarization effects have 



and higher-dimensional Randall-Sundrum models 

been studied for massive fermion fields in the context of the global monopole in |5fl. 

In this paper we will focus our attention on a new type of Casimir piston modeled on a conical manifold. 
The main difference of this configuration from the previous ones studied in the literature is the presence of 
a singularity at the origin and the fact that the piston is a curved manifold. As a result, only for particular 
dimensions the force will be unambiguous. This setup is particularly important in order to study the effects 
on the piston due to the presence of a geometric singularity. The main physical interest for the study of 



conical manifolds lies in field theoretical models requiring an orbifold compactification [28, 46]. In fact 
an orbifold is defined locally by the quotient space of a smooth manifold X and a discrete isometry group 
G. The action of the group on the manifold has, in general, fixed points which are mapped to conical 
singularities in the quotient space. This is a topic of great interest especially in the ambit of string theory 

I ID. 

We will utilize ^"-function techniques in order to obtain analytic expressions for the Casimir energy of the 
conical piston. The spectrum of a self-adjoint partial differential operator on compact manifolds is discrete, 
bounded from below and the eigenvalues A n form an increasing sequence. In this case one can define the 
spectral ^-function associated to the operator as 



(i.i) 



n=l 



which is convergent for %(s) > D/2, with D being the dimension of the manifold under consideration. One 
can analytically continue, in a unique way, £(s) to a meromorphic function with only simple poles in the 
whole complex plane which coincides with (11.11) in its domain of convergence. 

The outline of the paper is as follows. In section 2 we describe the geometry of the conical piston and 
we introduce the basic objects needed for our study. In particular, we consider two types of boundary con- 
ditions, namely Dirichlet and Neumann boundary conditions. In the framework of ^"-function regularization 



3 



we obtain expressions for the Casimir force which explicitly show their dependence on the particular ge- 
ometry of the piston. We specialize our general formulas to the case in which the piston is a <i-dimensional 
sphere and give very explicit results in particular dimensions. The Appendix contains a list of polynomials 
needed for the computation and the Conclusions point to the most important results of the article. 



II. THE CONICAL PISTON 



We will consider a particular manifold which is termed the bounded generalized cone. The generalized 
cone is defined as the D - (d + l)-dimensional manifold ^# = / x JV where jV is the base manifold, 
assumed to be a smooth Riemannian manifold possibly with boundary, and / = [0, 1] c R. The manifold 



is locally described by the hyperspherical metric jlol 



ds 2 = dr 2 + r 2 cTL 2 , (2.1) 

where dL 2 represents the metric on jV and r € I. It is known [7] that the curvatures on j$ and on the base 
jV are conformally related as follows 

R ij u = [R U ki - (d\5 j i ~ tfi&k)] , R'j = -id- l)<5'y] , (2.2) 

R = [R - did - 1)] , 

where R and R are the curvature tensors, respectively, on M and J/ ' . It can be readily understood, from the 
relations (12.2b . that, in general, the manifold under consideration has a singularity at the origin r = 0. 

For this type of singular Riemannian manifold the heat kernel and functional determinant of the associ- 



28] 



ated Laplace operator have been studied for massless and massive fields in ||7j, 

Let us next describe the piston configuration that arises from the generalized cone. Consider a cross 
section of the manifold ^# positioned at the point r = a with a e (0, b) and b > 0. We will denote this 
cross section by Jf a . Note that for any x and y in the interval (0, b), ^V x and jY y are diffeomorphic. The d- 
dimensional manifold Jf a naturally divides the manifold M into two distinct regions: region / represented 
by Mj = [0, a] X JV and region // represented by Mjj = (a, b] x jV where the first, Mj, contains the 
singularity at r - 0. Mj and Mjj represent two D-dimensional manifolds with boundary, where dMj = 
{0} U jV a and dM n = jY a U jV b . 

Clearly, the generalized cone M is obtained from the union of Mj and Mjj along their common boundary 
jV a , more precisely 
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The configuration that we have described above is a conical piston, where the piston itself is modeled by 
the cross section jV a of j$ at the point r - a. 

We will consider, in what follows, the Laplace operator on ^# acting on the space 5£ (./#) of 
square integrable scalar functions on the generalized cone. The starting point of our analysis is the following 
eigenvalue problem 

+ m 2 )(p = a 2 <p , (2.3) 

where we have introduced the spectral parameter (mass) m in order to make some of the subsequent integrals 

well defined. At the end of the calculation the limit m — > will be taken, giving results for massless scalar 

fields. In hyperspherical coordinates the Laplacian A ji takes the form of a Bessel type operator as follows 

d 2 d d ] 
or 1 r or r L 

with A,/)/ denoting the Laplace operator on the manifold jV . 

The idea is to solve the eigenvalue problem (12.31) in both region / and region //. The fields in one 

region are independent of the fields in the other region and therefore the corresponding spectral problems 

are independent. 

Region /, corresponding to the manifold Mj, contains the conical singularity at r = 0. We require the 
solution of (12.31 ) to be regular at the origin, and we obtain 

<pi = r^J v ( 7l r)®{n) , (2.5) 

where J v is the Bessel function of the first kind and we have set aj = y 2 + m 2 . The angular functions O(Q) 
represent the hyperspherical harmonics on JV satisfying the eigenvalue equation 

A = -i 2 <D(0) . (2.6) 

Region //, represented by the manifold Mjj, does not contain the conical singularity at r = 0. Therefore, 
a general solution to the eigenvalue problem (12.31 ) will be a linear combination of Bessel functions of the 
first and second kind as follows 

<p n - r¥ [A J v (y n r) + B F v ( 7// r)]0(Q) , (2.7) 

where A and B are arbitrary constants, and we have set a 2 n = y 2 u + m 2 . The index v of the Bessel function 
is easily found by substituting the general solutions (12.51 ) and (12.71 ) into (I2.3I ). By taking into account 
the angular relation (I2.6I ). one finds that the radial differential equation, in both regions, is satisfied if the 
following holds 

, o (1 - d) 2 
v 2 = A 2 + K - ' . (2.8) 
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We would like to introduce, at this point, the ^-function associated to the spectrum of the operator (12.41 ) 
in regions / and //. Since the form of the respective ^-functions is the same, we will utilize a unified notation 
and define 

m = 2>- + m 2 )~ s , (2.9) 

7i 

where i represents either / or // and we will assume that no negative eigenvalues occur so that we can use 
the standard branch cut of the logarithm. It is clear that the ^-function associated with the piston on the 
manifold ^# is obtained by adding the contributions from both regions, 

U(s) = + iu(s) ■ (2.10) 

In the following, we will express the ^-function on Mj and M]j in terms of the zeta function t^jy of the 
manifold jV \l6] defined as 

UU) = Yjd(y>~ 2s , (2.ii) 

V 

where d(v) is the degeneracy of the scalar harmonics <J>(n) on jV . This definition will serve the purpose 
of keeping the manifold JV unspecified throughout the calculations without jeopardizing the possibility to 
impose boundary conditions |7fl. 

It is useful to mention tha t in the framework of ^"-function regularization the Casimir energy is defined 

as follows HQ, HP 



23, 



34], 



£ Cas = • (2.12) 



a->0 2 \ 2 

where ji represents an arbitrary parameter with the dimension of a mass. By expanding the previous expres- 
sion about a - 0, one obtains a result which manifestly shows the structure of £cas, namely 

£ Cas - ^FPOf \-\ \ + \ I i + ln// 2 J Res U l-ij + 0{a) , (2.13) 

where Res denotes the residue of the function and FP its finite part. The last formula indicates that the 
Casimir energy has an ambiguity proportional to the heat kernel coefficient a^^ of the operator -A^ + m 2 , 
which, in turn, is given by the residue of ^ at s = -1/2. 

In the following calculations it will therefore be sufficient to compute the residue and finite part of the 
spectral ^-function on the manifold ^# at the point s = -1/2 to obtain £cas- When studying traditional 
pistons, only the finite part of the spectral ^-function depends on the position of the piston a. The force is 
then easily obtained from £c as by exploiting the formula 

FcM = ~^E Cas (a) , (2.14) 
oa 
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which gives an unambiguous result. The situation for the conical piston is, however, different. In this case 
the residue of the spectral ^-function at s = -1/2 will be, in general, dependent on the position a of the 
piston, as it will be shown later. The resulting force will be obtained from (12.131 ) in the form 

FcAa) = -\^ W ^\-\A ~ +ln ^) ^ Re&U (~h a ) + ° (a) ' (2 ' 15) 

It is thus clear that an unambiguous prediction of the force cannot be obtained in general in this setting. 
However, as we will see, if the piston jV is an even-dimensional manifold without boundary, the resulting 
force will be free of divergencies. 

In order to explicitly compute the spectral ^-functions in the two regions we need to impose specific 
boundary conditions, which, in turn, will provide implicit equations for the eigenvalues. In this work we 
will consider two types of boundary conditions, namely Dirichlet and Neumann boundary conditions. 



m. DIRICHLET BOUNDARY CONDITIONS 

In this section we will study the conical piston endowed with Dirichlet boundary conditions. In region / 
we impose Dirichlet boundary conditions on the piston Jf a positioned at r = a, which gives the following 
implicit equation for the eigenvalues yj 

Uyia) = . (3.1) 
In region // we impose Dirichlet boundary conditions on dMjj to obtain the linear system of equations 

A J v (y n a) + B Y v {y n a) = 
AJ v (y n b) + BY v (y n b) = . 

The above linear system has non-trivial solutions if the following equation is satisfied 

J v {yiia)Y v {y n b) - J v {y n b)Y v {y n a) - . (3.3) 



Since the index v is real and a is positive, equation (13.31 ) has only real and simple solutions, which implicitly 
determines the eigenvalues yu in region //. 

The starting point of our analysis is the representation of the spectral zeta function in terms of a 
contour integral in the complex plane valid for %(s) > (d + l)/2 |7-9, 25,(34]. By following the standard 
procedure, in region / we obtain 

1 C d 
Us, a) = d (y)^~ \ dk [ kl + m2 ]~ S dk ln [ rV ^ a) ] ' (3 " 4) 
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where m, as mentioned, is a mass which will be sent to zero at the end of the calculation, and T is a contour 
that encircles all the zeroes of J v (ka) on the positive real axis in the counterclockwise direction. The term 
k~ v has been introduced in order to avoid contributions coming from the origin in the subsequent contour 
deformation. Since the origin is not included in the contour the additional term will not change the result of 



the integral H34TI - In region II we have a similar representation, namely 

Oi(s, a) = ^ d(v)— dK [k 2 + m 2 \ S — In [J y (m)Y v (Kb) - J v (Kb)Y, 



(tea)] , 



(3.5) 



here F' is a contour that encircles all the zeroes of J v (Ka)Y v (Kb) - J v (Kb)Y v (Ka) on the positive real axis in 
the counterclockwise direction. 

By deforming the contour of integration in (13.41) to the imaginary axis and by performing the change of 
variables k — > vk/a we obtain for region I 



£i(s, a) = ^ d(v)d(s, a) , 



where 



sin(Tr^) r°° 
4 7 (s, a) = I dk 

JZ J ma 



2,2 



v l k 



m 



-ln[r v / v (v^)] 



(3.6) 



(3.7) 



An analogous deformation to the imaginary axis of the contour P in ( 13.51 ) and the change of variable k — > kv 
leads to the expression 



£n(s, a,b) = ^ d(v)C n (s, a, b) , 



(3.8) 



with 



C„(s, a, b) = f dK Uvk) 2 - m 2 ] 5 ^- In [K v (vKa)I v (vKb) - K v (vKb)I v (vKa)} . (3.9) 

71 Jm l 1 dK 

The results ( 13.71 ) and ( 13.91) are well defined in the strip 1/2 < 5l(s) < 1, and have been obtained by exploiting 
the formulas for Bessel functions of imaginary argument J3I L 

Let us for the moment focus our attention to the spectral ^-function in region /. In order to analytically 
continue the result to values of s for which %(s) < 1/2, we follow the methods of litis, 34] and utilize the 



asymptotic expansion of the modified Bessel functions I v (k) for v — > 00 and for z = k/v fixed. In detail 



[24, 



40] 



r v (yz) 



1 



0^1 



k=\ v 



&d+z 2 ) 1/4 

where the polynomials u^t) are determined by the recurrence relation 

1 



u k+ i(t) = -t 2 (l - t 2 )u' k (t) + 



If 

8 Jo 



dx{\ - 5r )uk(T) , 



(3.10) 



(3.11) 
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with uo(t) - 1 and 



t = 



r] = Vl + z 2 + In 



1 + Vl +z 2 



(3.12) 



The desired analytic continuation of 03.71 ) is obtained by adding and subtracting N leading terms of the 



uniform asymptotic expansion ( I3.10I ). One finds 17 



28, 



34] 



N 



fris, a) - Zi(s, a) + Aj 7) (s, a) , 



(3.13) 



i=-i 



where, once the limit m — > is performed, we have 



1 4 T(^+ 1) \ 2/ 

A ( ( , 7) (s,a) - ~U(s) , 

m a 2s ( r(j + & + |) 

Af\s, a ) = -— -U (s + = V *tf ^ 



(3.14) 
(3.15) 
(3.16) 



The remaining term, namely Zj(s, a), represents, by construction, an analytic function for %(s) > (d — 1 
AO /2 which is defined as 



J/Cs, ^ = J ^(v)Z^ 7) (5, a) 



(3.17) 



with 



7T J() d/c [ 



&d+^ 2 ) 1/4 

The terms D„(?) appearing in (13.131) are defined through the cumulant expansion [7- 



N 



In 



k=\ J n=l 



(3.19) 



and have the polynomial structure 



D n (t) = J ,/ !+2 ' 



(3.20) 



i=0 



For a list of the first polynomials D n (t) the reader is referred to the appendix. 

The spectral ^-function in region //, given by the integral ( 13.91 ), can be conveniently rewritten as a sum 
of three distinct terms 



C u (s, a, b) 



+ 



(tt£) r 

n Jm 

V 

Jl Jmb 



V [ V 2 « 2 2 

— : m 



cIk 
cIk 



V [ V 2 « 2 2 

— m 



b 2 



-In^fv*)] 



— In [k- v I v (vk)] + &Us, a, b) , 



(3.21) 
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where ^ZSs, a, ft) has the following integral representation, once the limit m — > has been taken, 



_ 2 sin(ns) f°° _ 2s d 

,a,b) = v dKK — In 

Jo 3k 



1 



K v (vi<b)Iy(vKa) 

K v (VKO)I v (VKb) 



(3.22) 



The first two integrals in (13.211 ) need to be analytically continued to values of s for which %(s) < 1/2. 
In complete analogy to the methods used for the analytic continuation of £i(s, a), we consider the uniform 
asymptotic expansion of the modified Bessel functions K v (k) for v — > oo and for z - k/v fixed, namely 



[24, 



40] 



K v (vz) 



-VI] 



1 + 



£<-i>' 



u k (t) 



k=\ 



(3.23) 



2v(l + Z 2)l/4 

By adding and subtracting N leading terms of the uniform asymptotic expansion (13.231 ) and (13.101 ). we 
obtain, for £n(s, a, b), the following expression 



N 



£ n (s, a, b) = Z n (s, a, b) + J*b(j, a, b) + ^ a| /7) (.s, a, ft) , 



i=-i 



where 



^b(j, a, ft) = J] d(v)^(j, a, ft) , 



(3.24) 



(3.25) 



Aj 7/) (s,a) = (-ly'Af^,^ + AfV,ft) , 
and Z//(5, a, ft) is an analytic function for %(s) > (d - 1 - AO/2 given by the expression 

Z H (s, a,b) - ^ d(v)Zy ir> (s, a, ft) , 

with 



(3.26) 



(3.27) 



Zy \s,a,b) = a v 



+ fc 2 Jy -2,sin(^) 



rf/r#f — Un [k K v {kv) - In J — - 

n Jo d/c ( [ V 2y (1 + a- 

!<££> r^^A/ln^/^l-ln 

7T Jo 5/f ( 



"7 
0)1/4 



n=l 
N 



D n {t) 



(1=1 



(3.28) 



4511 . As k approaches the 



J2irv(\ +k 2 ) 1 ' 4 . 

The function defined in (13.221 ) can actually be studied in more detail J26 . 
integral in ^^(s,a,b) represents a well defined function in the region %(s) < 1/2. As k — > 00, we utilize 
the uniform asymptotic expansion of the modified Bessel functions ( 13.10b and (13-231 ) to obtain 

Ky(vKb)I v (vKa) 



exp |-2v [r}(bK) - r)(aK)]} 



(3.29) 



K v {VKO)I v {VKb) 

Since a < ft, we have rj(aK) < r\Q}K) and the contributions of the expression (13.291 ) as k — » 00 are expo- 
nentially suppressed. These remarks allow us to conclude that J^,(s, a, ft) defines an analytic function for 
%(s) < 1/2 and we can safely set s = -1/2 in ^j){s,a, b). 
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In order to compute the Casimir energy for the configuration under consideration, we need to evaluate 
£i(s, a, b) and £n(s, a, b) in the neighborhood of s - -1/2. We would like to point out that for the explicit 
calculations that will follow it will be sufficient to subtract the first D terms of the asymptotic expansions 
(13.101 ) and (13.231 ) [34]. Therefore for the rest of the paper we will set N = D. In this way the functions 
Zj(s, a, b) and Zjj(s, a, b) in (13.181) and (13.281 ) are, by construction, analytic in the strip - 1 < ^R(s) < 1 /2 and, 
consequently, will not contribute to the residue of the pole of the ^-functions at s — -1/2. The evaluation 
of the functions A ( p(s, a) at the point s = -1/2 needs more care since the pole structure of £/ and is 
encoded in these terms. To obtain a systematic expansion, we set s = -1/2 + a and expand the resulting 



expressions for A^\s, a) about a = 0. From the general theory of spectral ^-functions Il30l l34ll we have the 
following expansion in terms of the variable a 



U(<* - 1) = U(-V + <V(-1) + 0(a 2 ) , = LAO) + at AO) + 0(a 2 ) , 



(3.30) 



41 1 a ~ 2~) = ^ Res u H) + fp ^ H 1 + ou> ' ■ 



and, for all d + 1 > i > 2, 



By utilizing the formulas above we are able to obtain 

- -!- [(in a 2 + 2 In 2 + l)^(-l) + £>(-!)] + 0(a) , 



A mi a A, a ) = -i^£tR-± 

■ 2 / a Ana Ana 



(3.31) 



(3.32) 



(3.33) 



A^la- -,a\ = — — ResZ>|--|- — 



1 



1 



n i 



i 



FP£jr\-- +lna z Res^ -- 



1 



+ 0(a) . 



(3.34) 



The term in (13.161 ) with i = 1 should be treated separately due to the fact that s = is a regular point of the 
spectral ^-function £jy(s). More explicitly one has 



A { ! ) \a--,a\ = - 



1 



1 



a 



1 



\6na 



U(0) 



i 



+ 



I6na 



For d + 1 > i > 2 one arrives, instead, at the following expression 



1 



1 



a>i I i - 1 

— Res^f 



2a sfn 



16 

£V(0) + £>(0) ( In a 2 + 2 In 2 - — 



i- 1 



+ O(a) . (3.35) 



1 



2a sfn 



+ ojj (in a 2 + y + 2 In 2 - 2) Res £ 



;-i\ (i-l 
,jY I — ^ — | + £2,Res 



+ 0(a) , (3.36) 



where y is the Euler-Mascheroni constant and we have defined the quantities 

0)1 = 1j Xi -p—( — TV 5 a > = L Xi >p—( — t\A p + ~T~ 



(3.37) 
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with *F(jc) being the logarithmic derivative of the Gamma function. From the definition (12.101 ). we have the 
following expression for C^( s ) f° r Dirichlet boundary conditions 



1 

a , a, i 

2 



+ 



Z/ I _ 2' a ) +ZlI \ 2' a ' b ) + ^ D \^' a '' 



ID/2] , s D , , 



(3.38) 



where [x] represents the integer part of x. 

According to (12. 13b . we need to extract the residue and finite part of the ^-function above in order to 
evaluate the Casimir energy. From equations ( |3.33b -( f3T36l > and, by recalling ( 13.261 ), we explicitly have 

Res^(-i«,*) = -i(i + ^) R es^(4) + ^|]^.Res^(^)-^(-l) 



1 1 

+ T7Z7^(°) + ^7Z^ Res ^ 
2 ynb 



z~ 1 



\6nb 



while for the finite part we get 



1 



1 



FP^(--, a,b\ = Z I \--,a) + Z II \--,a,b) + ^\--,a,b]- — 



2 

[D/2] 



1 



1 



1 



1 



(3.39) 



1 



FP £, r — + In a L Res ^ — 



+ 



+ 



+ 



Yj ^2,-FP U (— T^) + W2 ' ( ln « 2 + r + 2 In 2 - 2) Res U (^J^) 
OaReBO (^)] - ^ [(21n2 + l)^(-l) + £,(-1)] - lpp^> (-1) 



1 



^(0) + ^(0)(21n2-^ 



I6nb 

+ o)i (y + 2 In 2 - 2) Res 



i D 

— y 



2 -^nb 



/- 1 



+ f2,Res 



i=2 

i - 1 



i- 1 



(3.40) 



From formula (12.151) we can then write an explicit expression for the force on the piston when Dirichlet 
boundary conditions are imposed 



1 



1 



2 y[na 2 



2 

[D/2] 

E 

1=1 



1 



1 



1 



1 



1 



Aa 2 



(2 - lna 2 )Res£> (-;r|-FP£>(-- 



+ Q 2 iRes^> 



W2iFP^ 

2i - 1 



2i- 1 



1 I 2 
4^ 



- u) 2 i (4 - In a 2 - y - 2 In 2) Res 

[D/2] 



2i- 1 



Res I - ij - ^ w 2; Res 



2i- 1 



+ 0(a) , 
(3.41) 



where the prime indicates differentiation with respect to the variable a. The last term represents the ambi- 
guity present in the force in general. 
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IV. NEUMANN BOUNDARY CONDITIONS 



The calculational procedure to follow in order to compute £jg at s = -1/2 for Neumann boundary 
conditions closely resembles the one used in the previous section for the Dirichlet case. We will therefore 
describe only the few changes that are necessary (7H9[]. Imposing Neumann boundary conditions in region 



/ leads to the following implicit equation for the eigenvalues yj 

1 y^j Mayi) + ayiJ' v {aji) = , 
while for region // we obtain a linear system of equations for the unknowns A and B 

A [(^) Uayu) + ay n J v (ay n )] + B [( Y v {ay n ) + ay n Y' v {ay n )] = 
A J v {by u ) + bynJ'yibyn)] + B Y v (by n ) + by n Y' v {by n )] = 

This system of equations possesses non-trivial solutions if the following equation holds 



(4.1) 



(4.2) 



1 -d 



Jv(ayn) + ay n J' v {ayu) 



l-d 



Y v (byu) + bynY' v {by n ) 



l-d 



Y v {ayn) + ayuY' v (ayii) 



l-d 



J v (byn) + bynJ'y(byii) 







(4.3) 



47] 



For real v and positive a, the equation (14.11) possesses only real and simple zeroes if v > (d - l)/2a B18L 
and we restrict our attention to this case. 

In complete analogy with the Dirichlet case, we obtain the spectral ^"-functions in the regions / and //. To 
distinguish the results for Neumann boundary conditions from the ones for Dirichlet boundary conditions 
we use an upper index N. We find 



and 



tf(s,a,b) = Y J d{v)C\s,a,b) , 



(4.4) 



(4.5) 



with the integral representations valid for 1/2 < %(s) < 1, 



M sin(;rs) f°° 

C (s, a) = I dk 

JT Jim 



2 k 2 



m 



— ln[k- v (J3I v (vk) + vkI' v (vk))] , 
ok 



(4.6) 



and 



&\ s ,a,b) = r dK\{vK) 2 -m 2 Y I- In E v (K,a,b) 

" n Jm 1 1 OK 



(4.7) 



13 



In (14.71 ). for typographical convenience, we have defined fi = (1 - d)/2 and 



Ey(K, a, b) = \j31 v {aK) + aid' v {aK)\ \pK v {bic) + bi<K' v (bK)] - \j3K v {aK) + aKK' v (aK)] \/3I v (bK) + bKl' v {bK)} . (4.8) 

In order to perform the analytic continuation of the above integrals to the region ^R(s) < 1/2 for Neu- 
mann boundary conditions we need, in addition to the asymptotic expansion (13.101 ). the following one for 

i'(vz) HQ 



1 e vll (l +z 2 ) 1/4 



2/rv 



1 + 



k=l 



(4.9) 



where the polynomials Vk(t) are determined by the recurrence relation 



v k (t) = u k (t) + t{t l - 1) 



1 

-u k -i{t) + tu k _ x {t) 



(4.10) 



By exploiting (13.101 ) together with (14.9b . we obtain the following uniform asymptotic expansion, which will 
be useful in the subsequent calculations, 

~~ * M n (t,/3) 



In \fil v (kv) + kvl' v (kv)] ~ In | + k 2 ) 1/4 



V" 



(4.11) 



with M n (t,fi) defined by the cumulant expansion 

In 



34] 



M n (t,0) 



(4.12) 



k=\ \ k=l J] n=\ 

Moreover, the terms M n (t,0) have a structure analogous to the D n (t), namely 



M n (t,/3) = Y j Zi,n(j3)t 



n+2i 



(4.13) 



i=0 



where, in this case, the coefficients Zi, n depend on the variable j3 (see appendix). With the help of the above 
uniform asymptotic expansion we write ff(s) as a sum of the terms 

D 



tf(s,a) = W,(s,a) + £ A\ NJ \s,a) , 



(4.14) 



i=-i 



where one can prove that 



A ( ^ J \s, a) = A { !\{s, a) and A ( " J) (s, a) = -A^\s, a) , 



(4.15) 



furthermore for i > 1, once the coefficients Xjj, are replaced with Zi t b 171. la. 13411 . 



A^ J \s,a)=A ( P{s,a). 



(4.16) 
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The last term, i.e. Wi(s,a), is an analytic function for -1 < %(s) < 1/2 defined as 

Wi(s,a) = Yd(v)W y I) (s,d) , 



(4.17) 



with 



W v '\s,a) = a 2s v- 2s ^^- f dkk- 2s ^-{ln[k- v Q3I v {kv)+kvr v {kv))] 
t Jo ok { 



- In 



— rva + k 2 ) l/4 

2k 



-Z 



n=\ 



M n {t,P) 
v" 



(4.18) 



For the spectral ^-function g^' v (s, a, b) in region // we have a representation similar to the Dirichlet case 
(13 .21b . more explicitly 



& v (s,a,b) = 





poo 


dK 


\ V 2 K 2 


- m 2 


~ s d 




f ma 

V 




a 2 




dK 


sin^s) 


poo 


d.K 


r 2 2 




~ s d 




b 2 


- m 2 






J mb 


3k 



where ^\Xs, a, b) is represented by the following integral, once the limit m — > is taken, 



&Us,a,b) = v~ 2s 



sin(Trs) f 00 _ 2s d 

dK k — In A v (k, a, b) , 

n Jo ok 



with 



A v (K,a,b) = 1 



\fiI v (avK) + avKl v (avK)]\j3K v (bvK) + bvKK' v (bvK)\ 



(4.19) 



(4.20) 



(4.21) 



\J3K v {avK) + avKK v (avK)]\pi v (bvK) + bvJd' v (bvK)] 
The domain of analyticity of ^ v N {s,a,b) can be found by utilizing arguments analogous to the ones used 
for J^(tf, a, b) in the Dirichlet case. In fact, as k — > 0, the integral ( 14.201 ) is convergent for %(s) < 1 /2. As 
k — > oo, by using the uniform asymptotic expansions (14.1 II ) and (14.26b . we get 



1 - A v (K,a,b) ~ exp {-2v[rj(bK) - r/(aK)]} , 



(4.22) 



which is exponentially suppressed since rj{bK) > rj(aK) for a < b. Therefore, we can conclude that 
,^P v N {s,a,b) defines an analytic function for %(s) < 1/2. 

The analytic continuation of the integrals in ( 14.191 ) is performed exactly like in the Dirichlet case. In 
particular, we are able to write ^jf(s, a, b) in the form of a sum 

D 

Z%(s, a, b) = W n (s, a, b) + & N (s, a,b) + ^ ^'"^s, a, b) , (4.23) 



i=-l 



where 



(4.24) 
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and one can prove that 



A\ NJI \s,a,b) = (-\yA^'"(s,a)+A^' ,) (s,b) 



<NJ)t 



(4.25) 



The representation (14.231 ) is obtained by exploiting the uniform asymptotic expansion (14.1 II ) and the fol- 
lowing one 

" ~ ' , n M n (t,/3) 



y^(l^ 2 ) 1/L 



In [-pK v {vK) - vkK' v (vk)] ~ In 
which is obtained from the expansion for £y(w) in (13.231 ) and the one for K' v (vk) HQ 



(4.26) 



«=i 



2v k 

The function W//(5, a, ft) in (14.231 ) is analytic in the region -1 < *R(s) < 1/2 and has the form 

W n (s,a,b) = ^ d(v)wi rI) (s,a,b) , 

V 

where 

^— ^ f d K k~ 1s 4- 1 In \-BK v (vk) - vkK'Jvk)] - In 

t Jo <M 



(4.27) 



(4.28) 



im 9 c 9,. sin(^"i) 



^"'(I+k 2 ) 1 / 4 



in 



—e v \\ + k 1 ) 1 ' 4 
In 



-z 



(4.29) 



The evaluation of the residue and the finite part at s - - 1 /2 of the spectral ^-function for the conical 
piston endowed with Neumann boundary conditions proceeds along the same lines of the previous section 
on Dirichlet boundary conditions due, in particular, to the relations (14.151 ), ( 14.161 ) and ( 14.25 \ . For t*L in a 
neighborhood of s = -1/2, we obtain the following expression 

^(a- l -,a,b\ = W f I- 1 -, a) + W„ I -i a, b\ + & N W-, a, , 



+ 



22<' /) (4 + ^ a ) + Z A r ) H +a ^)- (4 - 30) 



1=0 v ' i=-\ 

The residue and finite part of the above spectral ^"-function are readily obtained leading to the result 

[D/2] 

/ M — 1 \ 

Anb 



R«#(-i «,*) = ^I + ^)Res^(-i) + ^ [ |\,Res^ 

1 ° ti- 



2i-l\ 1 

I- 1— — )-— 



(4.31) 
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and, for the finite part, 



= wA—a \ + Wu(-j,a,b \ + *J-±a,b j + 1 



FP^|-- + In a 2 Res ^ -- 



+ ~^ Z ^ 2 < FP ^ (^2^) + ^ ( ln a2 + ^ + 2 ln 2 - 2 ) Res ^> (~f^) 
+ fi*Res*> (^)j - JL [(2 In 2 + l)Cr(-l) + ^(-1)] + ^FP^ (-5 

- (I - ^ k- (0) + ^ (0) (21n2 - 2) 1 - tL LH0) + t 



+ (y + 2 In 2 - 2) Res £js 



i- 1 



+ Q,Res 



j- 1 



(4.32) 



where, in the above formulas, we have introduced the notation 



CO; 



\L Zi >p 1 — TV' n i = L^p—, — 7v X¥ \p + -r) • 



(4.33) 



which is identical to (13.371 ) once jc (jP is replaced by Zi, p . 

By exploiting the relation (12.151 ) we obtain the expression for the force on the piston when Neumann 
boundary conditions are imposed 



(2 - huriR.es, , (-Ij-FP^I-i 



2 V^fl' 2 



+ Q 2 ,-Res£> 



■ g w 2; FP U ~ ( 4 - ln « 2 - T - 2 In 2) Res (^2~~) 



2i- 1 



+ 4^(r ln ^ 2 



[D/2] 



/ n 2 ^ (2i-i 

Res + — 2^ & 2 ,-Res Cv \~~f~ 



+ 0(a) . 
(4.34) 



The expressions (13.411 ) and (14.341) represent the Casimir force on the piston, endowed, respectively, with 
Dirichlet and Neumann boundary conditions, which is valid in any dimension and for any compact and 
smooth manifold jV . It is important to stress that more explicit results can be obtained once the geometry 
of the manifold Jf has been specified. In addition, the analytic functions that appear in the results (13.411 ) 
and ( 14.341 ), namely Zj, Zu, Wj, Wu, J^£>, and J^jv can only be handled in a numerical way. Although the 
expressions for the force are written in terms of an infinite series in the eigenvalues v, in practice, one only 
sums finitely many terms in such a way to obtain a numerical result with a prescribed accuracy. 

These results also show explicitly the terms that are responsible for the ambiguity in the Casimir force. 
A closer look at these terms demonstrates that they are proportional to the heat kernel coefficients a^+i)/2~i 
of the manifold jV with < i < [(d + l)/2]. It is clear that this ambiguity in the prediction of the force 
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disappears if the manifold JY is even-dimensional without boundary. However, even if a suitable manifold 
is chosen for which the Casimir force on the piston is well defined, the formulas (13.411 ) and (14.341) do not 
allow to extract any information about its sign analytically, instead a numerical study is necessary; see 
Section VI. 

In the next section we will study the asymptotic behavior of the Casimir force on the piston when both 
a and b are large, and when a — > 0, i.e. the piston approaches the conical singularity. 



V. ASYMPTOTIC BEHAVIOR OF THE CASIMIR FORCE FOR BOTH a AND b LARGE AND FOR 

a — > 



In this section we will study the asymptotic behavior of the Casimir force for Dirichlet and Neumann 
boundary conditions when both the parameters a and b are large, namely b/a — > 1 . For this analysis, we 
will employ the arguments outlined in I45I for the case of concentric spheres. Moreover, we will describe 
in detail only the Dirichlet case, since the Neumann case can be treated in exactly the same way. In the 
limit b/a — > 1 the Casimir force for the conical piston will reproduce the one for parallel plates. From 
the expressions (13.411 ) and (14.341 ) it is not difficult to realize that all the terms but and j^yy vanish 
proportionally to a" 2 as a — > 00. Therefore, the asymptotic behavior for a and b large when Dirichlet or 
Neumann boundary conditions are imposed is given, respectively, by and j^yv- With the last remark in 
mind we can write that, when b/a — > 1, 

K v (vKb)I v (vKa) 



1 c 00 d 

"c D >.»~-s2>»'J 



1 



(5.1) 



K v (vKa)I v (vKb) 

The derivative in the integrand, denoted by P(k, a, b), can be evaluated my making use of the properties of 



the derivative of the modified Bessel functions to get 

K v (vtcb) 



P(k, a, b) = 



aKy{vKa)I v {vKb) 



1 - 



Ky{vKb)I v {vKa) 
K v (vKa)I v (vKb) 



-1 



(5.2) 



By performing the following change of variables k —> K/a, we obtain the expression 

The b — > a behavior of the Casimir force is already captured using the uniform asymptotic expansion of the 
modified Bessel functions to obtain 14511 

Piit, a, b, n) 



(-. \ 00 00 

V A 1, *)~ 2v vT77£ e-H^ h i)-M J] 
a a ' n=l i=0 



V 



(5.4) 



18 



where pi are polynomials in t which vanish when b — > a and such that po = 1 14511 . By utilizing the inverse 
Mellin transform we can rewrite the sum of the exponentials in (15.41) in terms of a complex integral as 
follows 



(5.5) 



£ ^v^)-**)] = ^ _L f Ja r(a)(2nv)- 

which is valid for %(y) > 0, and we assume that 21(c) is large enough so that the sum over v and the integral 
in a can be safely interchanged. 

By substituting the expressions (15.41) and (15.51) into the integral (15.3I ). and by recalling the definition of 
gjy in (12.111 ) we obtain 

F^ia, b) ~ ^2 2 J ^r(a)(2)-^ R (a)^ (— - l) 



i=0 



/-•OO 

J dK^Jl + K 2 
JO 



1\* 



Pi(t, a, b, n) , 



(5.6) 



where the integral over k is convergent for 51(c) > 2, since T](bic/a) - tj(k) — > ln(ft/a) as k — > and 



T}(jbK/a) - t]{k) — > (ft /a - l)/c as a- — > oo [45]. 

It is convenient, at this point, to introduce the variable g = ft/a - 1. In this new variable the limit we are 
interested in corresponds to q — > 0. The leading behavior in a of the integral (15.61 ) is given by the term with 



i = 0, since the other polynomials pi, with i > 1, vanish as a — > 0, more explicitly 

~ ^2~2 J dar(ar)(2)-° , &(ar)^> (- - lj J <k Vl + * 2 ((<7 + 0*) - 'K«)]"° ■ (5-7) 

By closing the contour to the right, we encounter all poles of ^j/ with %(a) > 2. The rightmost pole of 
t,jr results from a = d + 2, the other poles further to the left will be irrelevant for what follows because the 
leading behavior of (15.71 ) as q — » is proportional to the residue of gjs(a/2 - 1) at a = d + 2. This is seen 
by using [45] 



j] ((<? +\)k)- n(K) = q Vl+* 2 + 0(q 2 ) . (5.8) 

The last remark, and the expansion ( 15.81 ), allow us to obtain the Casimir force for the conical piston 
endowed with Dirichlet boundary conditions when b/a — > 1 as follows 



4-' 



(5.9) 



where the higher orders in q would involve the other heat kernel coefficients of J/ '. The Neumann case is 
treated in exactly the same way and the resulting Casimir force when b/a — > 1 reduces to ( 15.91 ). 
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Let us now turn our attention to the limiting behavior of the Casimir force when the piston approaches 
the conical singularity, namely, when a — > 0. From the expressions (13.411) and (14.341) it is not difficult to 
observe that all the terms, except and J^jv> are proportional to cT 2 . The proportionality coefficient 
and its sign will depend on the specific geometry of the piston JV . It is therefore sufficient to study the 
asymptotics of J^zj and J?jv as a — > 0. Once again, we will focus on the Dirichlet case since the Neumann 
case can be treated in exactly the same way. 

From the expressions (13.221 ) and (15.21) . by exploiting the change of variables k — > k/v, we can write 

~\^(~h a,b )^~h^ d{y)y L dKp (^ a ' b ) • < 5 - 10 ) 

Since we are interested in the small a expansion, we utilize the series representation of the modified Bessel 
functions to obtain the following expression, to leading order in a, of (15.21 ) 

:) 2V + 0(k 2 - 2 ) . (5.11) 



T 2 (y) I v {Kb) \2j 



This expansion allows us to write 



1 \ 2 ^ a 2v - 1 r°° K v (Kb) (k\ 2v i 7v+0 , 



The integral in the last expression is convergent, since 

KyQcb) Ik\ 2v 1 
I v (Kb) \2 

and 



(^) ~ -T(v)T(v + \)b~ 2v , for ^^0, (5.13) 



K v (Kb)(K^ m _ 2 for (514) 



I v (Kb) Ml \2, 

We can therefore conclude that for v > -1/2, which is within the assumptions of our work, the contribution 
is subleading as a — > 0. 

The previous arguments show that as a — » the Casimir force on the piston possesses a behavior of the 
type aT 2 with the proportionality coefficient depending on £jy and, hence, on the geometry of the piston. 

In the next section we will study a particular case in which the manifold JV is a J-dimensional sphere. 
This case is of special interest because the spectral zeta function of a ^-dimensional sphere can be explicitly 
evaluated in terms of a linear combination of Hurwitz ^"-functions O, 



15, 



34]. 



VI. A SPECIFIC PISTON: THE cZ-DIMENSIONAL SPHERE 



In this section we will study the case in which the base manifold is a J-dimensional sphere. In this 
particular situation the eigenvalues of the Laplacian on JV are known to be 
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with / > 0, and the eigenfunctions are hyperspherical harmonics with degeneracy 

* ) = (2 ' +rf - 1) i|^if • (62) 

The explicit knowledge of the eigenvalues v and their degeneracy allows us to write the ^-function on Jf 
as follows 

-2s 



t (/ + rf-2) ! /, d-1 

The previous expression can be conveniently rewritten in terms of a linear combination of Hurwitz £- 
functions, namely |7, 80 



/ J — 1 \ 

= 2 J] e a & 2j - or - 1, — , (6.4) 

a=0 \ Z / 



where the coefficients e a are defined through the relation 

(/ + </- 2)! /, d-1 



/v=n \ 



/ + ^— . (6.5) 



As one can easily see, from the general results (13.411 ) and (14.34b . we need to compute either the value or 
the residue or the finite part of the spectral ^-function on the manifold JV^ at sp ecific points. For s = —ml 2, 
with m > -1, by utilizing the equation (16.41) . we obtain the result |2, 

d-\ 



a=0 

where we have used the relation, valid forn > |3I ], 



/V — I I x ' 



y , , B n+ i(q) 

£ H (-n, q) = — , (6.7) 

n + 1 



with B n (q) being the Bernoulli polynomials. It is clear from the relation (16.61 ) that the residue of £jy(s) at 
s = -1/2 and s - 1/2 vanish. 

The remaining values of s that we need to analyze are the ones of the form s = m/2 with m > 0. It is 
well known that the Hurwitz ^-function, g), has a simple pole at 5 = 1. The explicit form of its Laurent 
expansion depends on whether q is an integer or a half-integer. In fact, for q = n + 1/2, with n > 0, we have 



( H ll +sx,n + ij = -L +y + 21n2 - 2^ ^-y + O(x) , 



(6.8) 



while for g = n and n > we obtain 

1 1 

6Kl + «c,n) = — + y-Y- + <9(x). (6.9) 
sx f—fk 

k=l 
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By utilizing the above results and the relation (16.41 ) we obtain the following expression for the residue of 
£js(s) at s = m/2 with d > m > 2, 



in 



Res CyV [ -z ) = e m -2 , 



(6.10) 



the contribution to the residue coming from the index a = m - 2 of the summation in (16.4I ). According to 
(16.8b and (16.91) . for the finite part at s = m/2 we have to distinguish between two cases. If the dimension d 
of the manifold ,jY is even, then 



d-l i 

FPCr (f ) = 2 2 lm-a-1, 
while if J is ocW, we have 

FP^(y) = 2 ^ ««& 



j - r 



+ 2e m - 2 y + 2 In 2 - 2 V — 



(6.11) 



a=0 
at-m—7. 



m — a — 1 , 



) + 2 ^" 2 r "?2l^T 



(6.12) 



The results that we have obtained in this section so far are well suited for the analysis of conical pistons 
with Dirichlet boundary conditions. However, the case of Neumann boundary conditions needs a special 
treatment when the manifold jY is a J-dimensional sphere. In fact, from the integral representation (14.61 ) of 
the spectral ^"-function in region /, one can notice that the case 1 = 0, namely v = (d— l)/2, needs particular 
care because the behavior of ^/3Itti(ak) + akl' d _ x (ak)j as k — > is different. This corresponds, according to 
(16.11) . to the lowest angular eigenvalue v. In this case, in fact, we have the following small k expansion 

(ak) 



[pLj^iak) + akl'^(ak)j = ^ ' +0[k^ 

2 2 r {~) 



(6.13) 



The above expansion suggests that a more suitable integral representation for the contribution of the mode 
v = (d-l)/2is the following 



m i-n t c / d - 1 \ sin(7rs) r°° 9 . d \ d+3 I 

tf' l -°(s,a) = a 2i \—- 1 i - J dkk- 2i — ln^k-—[/3Id^(k) + kI' d _ i (k) 



(6.14) 



which is valid in the strip 1/2 < %(s) < 1. In order to analytically continue the above expression in the 
neighborhood of s = - 1 /2, we utilize the following k — > oo asymptotic expansion 



In 



V 2 2 



In 



V2^ 



E 

n=l 



it" 



(6.15) 



where the polynomials are defined according to the relation 
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The expressions (16.151 ) and (16.161 ) can be obtained by exploiting the asymptotic expansions for z — > °° [39], 



1 + V^ 



and I'Jz) 



V2^z 



i + y Pn{y ^ 

Zj z n 

n=\ 



(6.17) 



where the functions r n (x) and p n (x) axe obtained from the relations 



l r 2 i 1 42,421 n 

r„(.v)=l. n(,) = =^, r„(x) = 4 * tg f] [4^ 2 - (2/ - 3) 2 ] , (6.18) 

1=2 



«!8" 



and 



1 " 

- 1 . p«w = ^ n [ 4x2 ~ (2/ ~ i)2 i 

i=\ 



(6.19) 



At this point it will be sufficient, for our purposes, to add and subtract the leading term of the expansion to 
obtain the result which is valid for -1 < %(s) < 1/2, 



tf J -°(s,a) = W' [ - u {s,a)-a 



1=0, 



2s sin(7rs) Id — l\ ( Id \ 1 1 



2s 2s — 1 



1 (3 d-\ (d-Y) 1 
- + — + — 



2s + 1 \8 



(6.20) 



where 



W'f\s,a) = a 2s (^)^^^^r 2 ^|ln[r^^/ ¥ W + H^w) 



H(k - 1) 



In 







dkk~ z 


-< 


dk\ 


All 


d 




+ — 







(6.21) 



with the step function H(k). In region // a similar argument for the lowest angular eigenvalue is necessary. 
In fact by utilizing the small k expansion of the modified Bessel functions, it is not very difficult to show 
that Ed-i (a, k) in (14.81 ) behaves differently as k — > 0. More specifically, one obtains 



2a 2 



d+3 d-l 

a\ ~t~ (a\~ 



(6.22) 



A suitable representation for the spectral ^-function in region // corresponding to the lowest eigenvalue is 
therefore 



,nj=o, s 2s (d- \\sin(jrs) 
Zj/ (s,a) = a 



+ 



+ 



2 S (d - 1 \ sin(ns) 



n Jo ok 

^1 r *,-#i„ 

7T J,, OK 



K^i-pK^-KK'^K)) 

2 2 

d+3 , 
K- — {J3h_ 1 {K) + Kl' d _ l {K)) 



(d-l\ S m(ns) C°° _ 2s d 

dKK — In Arf-i (K,a,b) . 

\ 2 / n Jo 6k 2 



(6.23) 
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After analytic continuation of (16.231) . we obtain the following result which is well defined in the region 
-1 < %(s) < 1/2 



^ J -°(s,a,b) = W'jj\s,a,b) + &'„ u (s,a,b) - 



rl=0 



sin(Ks) Id — l\ ( 1 



IX 



2s 



1 



+ 



2s + 1 



3 + d-l + (d - \y 



a ls - b 2s 
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where Wjj ^, a, b) has the form 
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Let us mention that ^j^°(s, a, b) coincides with (14.201 ) once we set v - (d - l)/2. 

From the above results, it is not difficult to obtain the residue and finite part for qjz '~° at s = - 1 /2 
corresponding to the lowest angular eigenvalue. In more detail, we have 
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Therefore, the contribution of v = (d - l)/2 to the Casimir force on the piston when Neumann boundary 
conditions are imposed is, according to (12.15I ). 
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Obviously, the total Casimir force on the piston for Neumann boundary conditions is the sum of (16.281 ) 
and (14.341 ) where in the latter the lowest angular eigenvalue is omitted. In the next subsections we present 
explicit results for specific dimensions d. We would like to point out that in the formulas that will follow it 
is understood that the functions Z, W and & are evaluated for the specific dimension under consideration. 

A. Specific Dimensions for Dirichlet Boundary Conditions 



In the following special cases, we will set, for simplicity, 6 = 1. When the piston jV is a sphere of 
dimension d = 2, and, therefore, the dimension of ^ is D - 3, we obtain 
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For d - 3, or D - 4, we have the result 
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For d = 4, or D = 5, we obtain 
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And, finally, for J = 5, or D = 6, we get 
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As mentioned previously we can see, from the above results that the Casimir force on the piston is not a 
well defined quantity when the piston jV is odd dimensional. The Casimir force on the piston for d = 2 and 
d = 4, is shown in Figure [TJ We can see that, for d = 2, the piston is repelled from the conical singularity 
and attracted to the base manifold positioned at b = 1. For d = 4, instead, a point of unstable equilibrium 
is present. If the piston is to the left of this point, it is attracted to the conical singularity. If it is to its right, 
then the piston is attracted to the base manifold. 



(a) d = 2 (D = 3) 



(b) d = 4 (D = 5) 



Figure 1: Plots of the Casimir force, F®" s {a), on the piston jV for Dirichlet boundary conditions as a function of the 
position a. 



B. Specific Dimensions for Neumann Boundary Conditions 



For Neumann boundary conditions, we need to consider the results obtained in (14.341) . with the lowest 
angular eigenvalue omitted, and the one obtained in (16.281) . It is useful to stress that the spectral ^-function 
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Cyy(s) with the lowest angular eigenvalue omitted is related to the ordinary £yf (*) through the relation 

I d — \ \ -2s 

u(s) = uw- HH ■ ( 6 - 33 ) 



This expression is utilized in the results that will follow in order to obtain the Neumann Casimir force. 
For d -2, thus D - 3, we have the following expression for the force 
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For d = 4, thus D = 5, we obtain 
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Once again, we would like to point out that also for Neumann boundary conditions the Casimir force 
on jY is not well defined when d is odd. The Casimir force on the piston for d = 2 and d = 4, is shown 
in Figure |2] We can notice that in both cases, namely d - 2 and d - 4, there exists a point of unstable 
equilibrium. If the piston is to the left of this point it is attracted to the conical singularity, while, if it is on 
its right, it is attracted to the base manifold. 
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(a) d = 2, and D = 3 (b) d = 4, and £> = 5 

Figure 2: Plots of the Casimir force, F^(a), on the piston Jif for Neuman boundary conditions as function of the 
position a. 

VII. CONCLUDING REMARKS 

In this paper we have analyzed the Casimir effect for massless scalar fields in the setting of the bounded 
generalized cone in the cases in which the field is endowed with either Dirichlet or Neuman boundary 
conditions. The piston geometry has been obtained by dividing the generalized cone into two regions 
separated by its cross section positioned at a with a e (0, 1). The cross section itself represents the piston, 
and the structure so obtained has been denoted conical piston. By utilizing ^"-function regularization meth- 
ods, we have obtained very general results for the Casimir energy and, hence, for the corresponding Casimir 
force. The expressions obtained for the Casimir force are valid for any piston jV and any dimension D and 
are written in terms of the spectral ^-function In order to obtain more specific results, we have consid- 
ered in detail the case when the manifold JV is a J-dimensional sphere. Explicit numerical results for this 
case have also been obtained for d = 2 and d = 4. 

We would like to stress that this work provides a study of a piston geometry that is more general than 
the ones considered in the literature where the two chambers have fundamentally the same type of geom- 
etry. The two separate regions of the conical piston, instead, have different geometries since one of them 
contains a singularity while the other does not. This is the feature that makes the conical piston particularly 
interesting. These types of conical pistons are important, e.g., for the study of the interaction, through the 
Casimir force, of a piston and a singular point (such as the tip of the cone). This results may also have some 
relevance in in the framework of field theories with orbifold compactification. 

The next step in this kind of investigations would be to analyze the Casimir force for massless scalar 
fields on the conical piston satisfying mixed, or hybrid, boundary conditions. It would be interesting to 
further understand how the Casimir force between the piston and the singular point of the manifold 
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is influenced by the boundary conditions imposed. Furthermore, it would be desirable to have a deeper 
understanding of the influence of the geometry on on the force. 
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Appendix A: Polynomials D„(t) and M n {t,j3) up to the order n = 6 



the sixth order. 



In this appendix we list for completeness the polynomials D n (t) and M n (t,f5) up to the sixt 
By utilizing the cumulant expansion ( 13.191 ) and the recurrence relation (13.111) one obtains 
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The polynomials M n (t,fl) follow by using ( 14.121 ) together with the recurrence relation ( 14.101) . One can 
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